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I. INTRODUCTION 

One of the most outstanding difficulties of experimen- 
tal light quark spectroscopy - like in studies of char- 
monium decays to light quark mesons at BES III [T] 
or future studies of photoproduction at GlueX - is in 
the disentanglement of overlapping and interfering me- 
son states, which often have widths of several hundreds 
of MeV. This requires amplitude analyses, where exper- 
imental distributions are described by a seies of theo- 
retical amplitudes ( decay amplitudes ) with each am- 
phtude generally multiplied by a freely fit parameter ( 
production amplitudes). In the past, decay amphtudes 
were generally written using the isobar model, i.e. as- 
suming a multi-particle decay proceeded through a se- 
ries of two-body resonance decays with the resonance de- 
cays usually parametrized as Breit-Wigner amplitudes. 
This model, however, is known to violate unitarity. With 
high-statistics data samples now available at BES III and 
later in GlueX, as well as other current and future exper- 
iments, more careful attention must now be paid to the 
theoretical descriptions of the decay amplitudes, and phe- 
nomena such as final-state re-scattering and inelasticity 
must be considered. 

The decay J/ip — > Tr+Tr^Tr", which is observed to pro- 
ceed dominantly through pn, provides a simple context 
in which re-scattering effects can be studied. Here the 
TTTT system is limited to either J^^ — 1 (P-wave) or 
3 (P-wave). Neglecting the small 3 component, this 
reaction thus provides clean access to P-wave tttt scat- 
tering. The decay J/ip — tt'^tt^tt'^ has previously been 
studied experimentally by BES II [5] and BaBar 3 , but 
limited statistics prevented any detailed analysis of the 
3tt substructure. BES III will soon have a set of J/ip de- 
cays many times larger than what is now available, and 
this data set could be used to greatly improve many of 
the theoretical uncertainties associated with re-scattering 
effects. 

In this work, we present a coupled channel analysis of 
J/ip — >■ 7r"'"7r~7r'' decays in which we consider both tttt 
and KK isospin-1 intermediate states. In particular, we 
take advantage of unitarity constraints to reconstruct the 



amplitudes based on their analytical properties. Unitar- 
ity relates the discontinuity of the isobar amplitude to 
the scattering amplitude and we use the available data 
on P-wave tttt scattering to construct analytical tttt and 
KK scattering amplitudes. We show that available data 
on the 37r decay of the J/ijj is inconsistent with the single 
channel parametrization. The effect of the intermediate 
KK pairs is to enhance the contribution from the tail of 
the p(770) while reducing contributions from higher-mass 
p excitations. 

This paper is organized as follows. In the following 
section, we discuss the analytical properties of the pro- 
duction and scattering amplitudes. We also construct 
an analytical model for two-channel tttt and KK scatter- 
ing and finally compare theoretical predictions with the 
experimental data. A summary is given in Section [lll[ 

II. P-WAVE TTTT EFFECTS IN J/iP -> Tr+Tr^Tr" 
DECAY 

For each helicity state, A, of the J/ip, the amplitude to 
decay to three pions is a function of three angles and two 
invariant masses. In the rest frame of the J/ip, the angles 
may be chosen to specify the orientation of the plane 
formed by the momenta of the three produced pions with 
respect to the direction of polarization of the J/ip. The 
invariant masses correspond then to the Dalitz variables 
describing the 37r system. Denoting the four-momenta 
by p±.o, P for tt^, and J/ip, respectively, the general 
expression for the amplitude is given by 

{TT\+7T-,out\J/i;{\),tn) - {2nfSH ^ ft - P)tT^, 

1=0, ± 

(1) 

with, in the rest frame of the J/ip, 

Tx = -ie{X) ■ (p+ X p_)P(s+_, S0+, s-o). (2) 

Here e is the polarization vector of the J/ip, the Dalitz 
invariants are defined by Sij = {pi ~\- pjY for i,j = ±,0 
and satisfy s+_ + so+ + s-o = M"^ + 3ml, pi = Pi/|pi|, 
and the scalar form factor P describes the dynamics of 
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the decay. It is |Fp that determines the distribution of 
events in the Dahtz plot, i.e. |Fp = const, yields a flat 
distribution. Since Pi = in the J /if) rest frame, any 
two pion momenta can be used instead of p-|- and p_ in 
Eq.([2]) to specify the orientation of the decay plane. 

The isobar model makes a specific assumption about T, 
i.e. the decay is assumed to proceed via a quasi two-body 
process in which a pair of pions in a low partial wave and 
a spectator are formed without any further interactions. 
The isobar model violates unitarity, which forces interac- 
tions between pions from the quasi two-body state and 
the spectator to be included. If the quasi two-body state, 
however, is dominated by a low-mass, narrow resonance, 
then the overlap between the resonance and the specta- 
tor pion wave functions is expected to be small. Indeed, 
in the case of the ttttN final state at a total center of 
mass energy below 2 GeV [H [5] (one of the very few phe- 
nomenological analyses of re-scattering effects in three- 
particle systems that we are aware of), the re-scattering 
corrections were found to not exceed 20% [S| . In the case 
of the J /ip with even higher center of mass energy and 
with a pronounced p resonance in tttt, we expect these 
effects to be even smaller. Nevertheless, it will be impor- 
tant to quantify the size of such re-scattering effects in 
three-body J /if) decays, in particular in view of the very 
high statistics data currently being collected at BES III. 




q^+q_=0 



FIG. 1: Definition of the decay angles in the J/ip rest frame 
(left and top right) and the ■k'^-k~ isobar rest frame (bottom 
right). 

The two lowest tttt partial waves allowed in J /ij} decay 
have L = 1 (P) and L — (F). Little is known about 
higher partial waves, but the F-wave is already very weak 
with the phase shift staying below 5" for energies up to 



1.45 GeV fT". In the following we will thus keep only the 
P-wave in our isobar analysis. Within the isobar model 
with a single P-wave tttt isobar, the amplitude T in Eq.([T]) 
is given by 

= E E Dl:^i^Ml^M)Ms3k) (3) 

i=0,± M=±,0 

where the angles are illustrated in Figjl] and the indices 
ijk run through cyclic permutations of 0, — 0[9]. Here 
A is the spin projection of the J/'0, which, together with 
the X and y defined with respect to a lab coordinate sys- 
tem, defines the z axis. The rotation is given by three 
Euler angles, = rk{4>k, i^k^ipk), which rotates the stan- 
dard configuration that corresponds to the {ij)k coupling 
scheme (with the ij forming the L = 1 isobar and tt*^ be- 
ing the spectator) to the actual one. In the standard 
configuration tt*"' has momentum along —z and tt' and tt^ 
have momenta in the xz plane with tt' having a positive 
X component. Finally, 6k is the polar angle of the tt* 
in the tt^tt-' rest frame. In other words, (j)k and are 
the azimuthal and polar angles, respectively, of the total 
momentum of the ttV^ pair in the Stt rest frame, while 
V^fc and Ok are the azimuthal and polar angles, respec- 
tively, of the TT* in the ttV-' rest frame {i.e. the isobar 
rest frame). For the three possible coupling schemes, the 
corresponding Euler rotations, r^, z = ±, 0, are related to 
each other by 

ro =r+r(0,x+,0) =r_r-i(0,X-,0), (4) 

where x+(X-) is the angle between tt"*" (tt^) and tt" in 
the Stt rest frame. This enables us to write T in terms of 
D{ro) alone: 

+ <.(x+)4.o(^+)^m(s-o) + <,(X-)<,o(^-)J^m(so+)] . 

(5) 

The helicity amplitudes, F^, are linear combinations of 
the L — S coupling, isospin-/ amphtudes, P/^g [10 . In 
the case considered here with I ^ L = S = 1, only a 
single amplitude, Px^i' contributes, and 

F^is.j) ^ ^^^{l^i\l,^i■,l,o)Fl,,{s,,), (6) 

which implies Fq = and Fi = — P_i. Finally, compar- 
ing with Eq.([2]), in the isobar model we obtain 



P(s+_, so+- s-o) 
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E('^-.i + '5.,-i)dLo(^o)i^i(.s+-) + (rft.(X-) + dl,,{x-))dloie-)Fi{so+) 



(7) 
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A. Unitarity constraints on the isobar amplitudes inelasticity, 77, these amplitudes are given by, 



Writing the J/V' decay amplitude as an analytical func- 
tion of the channel sub-energy, Sjk, one finds 



776 



1 



rje 



1 



2ipK 



^ttK ^Ktt 



(11) 



{{ij)k,out\J/i;,in) - {{ij)k,in\J/TP,in) ^ {2Tryi X ^^/j:^ 
'x^S^iPi+Pj -p[-p'j)t*{ij;i'j'){ii'j')k,out\J/ip,m), ^here the phase space factors are given by Pa{s) = 

* •' ^1 — Sa/s and the pa in Eq.M are defined as Pq(s) = 

(8) 4:q'^{s) Pais) — (s — Sa)Pa_is). Similarly one finds 



where t{ij;i'j') is the scattering amplitude between the 
incoming \ij,in) and the outgoing \i'j\out) state. The 
two matrix elements on the l.h.s. give the J/ip decay am- 
plitude evaluated at Sij+ie and Sij—ie, respectively. Sim- 
ilarly, discontinuities across the other two sub-channel 
energies can be considered. However, because of the sym- 
metry of the isobar amplitude under permutation of the 
three pious, they all lead to the same unitarity relation. 
The summation over intermediate states on the r.h.s. 
should include inelastic channels. It is known that the P- 
wave TTTT amplitude is elastic up to energies ^1.4 GeV, 
with the KK channel effectively saturating inelasticity 
above this energy, at least up to ^ 1.9 GeV where data 
is available. Thus, using a single KK intermediate chan- 
nel, Eq.([8|) leads to 



ImFjr{s 



il^(s)p,{s)F^{s)e{s-^ml) 
ilK{s)PK{s)FK{s)6{s - 4mj^). 



(9) 



As discussed in Section [TTj this is an approximate rela- 
tion, which ignores contributions to the r.h.s. from re- 
scattering between a pion from the isobar and the spec- 
tator pion. In Eq.(|9|, the helicity-1 isobar amplitude, 
Fi from the r.h.s. of Eq.Q, is denoted by Fjr{s) to dis- 
tinguish it from the corresponding helicity-1 amplitude 
for production of KK P-wave pair in J/-0 — {KK)pTr, 
which we denote by Fk{s). Furthermore we define Fa{s) 
(a — t:,K) as the reduced isobar amplitude, i.e. the 
amplitude with the angular momentum barrier factors 

2qa{s) = sjs - Sa, Sa = 4to^ , 

2p(s) = 



M2 



(10) 



removed, so that F^ = Fa/{2qa2p). Here Qq, is the rel- 
ative momentum between the pious (a = tt) or kaons 
(a = K) in the isobar rest frame, and p is the break-up 
momentum of the J/^p (mass M) into an isobar of mass 
^/s and the spectator pion. In addition, t^T^ (i-KR) is 
the elastic, isospin-1 tttt [KK] P-wave amplitude, and 
It^k is the P-wave transition amplitude for KK — tttt. 
Similarly, tap are defined as the scattering amplitudes 
without the barrier factors, i.e. taf} = tap I ('^<laq_f})- In 
terms of the P-wave phase shifts, and 5k, and the 



ImFK{s + ie) = i*K,{s)p^{s)F^{s)e{s-Aml) 

+ i*KK[s)pK{s)FK{s)0{s^Am\). 

(12) 

In the isobar approximation the form factors P^- and F^ 
are real analytical functions (Pq(s*) — F*{s)) of a single 
sub-channel energy and thus have only the unitary cuts 
and satisfy 



Pa(s) 



ImFa{s') 



ds'. 



(13) 



With ImFa given by Eqs.(^9^ and ( 12 1 the isobar form fac 



tors become a set of two coupled integral equations. An 
analytical solution can be obtained using the standard 
Omnes-Muskhelishvili approach [11] [T^] . To this extent 
one first notices that, in the two-channel {a = tt, K) ap- 
proximation considered here, the unitarity condition for 
the reduced scattering amplitudes, tap, is given by 

Imiapis + ie) ^ ^ i*a^{s)pj{s)9{s - s^)^l3{s). (14) 

This implies that the right hand discontinuity relations 
for Fa are satisfied by the functions [13] 



(15) 



where the production amplitudes, Pq.(s), are real for 
s > and free from right hand side discontinuities. If 
Pq(s) is to be free from discontinuities for s < then the 
production amplitudes Pq(s) have to satisfy the integral 
equation 



1 

TT J-c 



,ImPa{s') 



ds' , 

is' - s) 



(16) 



For s < 0, IniPa{s) is obtained from the condition 
ImFais) = 0, 

ImPa{s)= [Rei{s)]-^p[Imi{s)]p^ReP^{s). (17) 



In general, at most one subtraction in Eq.(16) may be 



needed based on the asymptotic behavior of the scatter- 
ing amplitude, which is discussed below. The subtraction 
constants would then become fit parameters in this uni- 
tarized isobar approach. 
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P-wave nn scattering amplitude: general 
properties 



In order to solve Eqs.(13l and (161, it is convenient 
to separate the left (s < 0) and right (s > s^) cut con- 
tributions to the reduced scattering amplitudes <q^(s). 
This can be done using the "N/D" representation inde- 
pendently for the amplitude of each channel [Hj , 



tap 



(18) 



with Nai3 — N^a and Dap — Dpa having only the left 
and right hand cuts, respectively. Then analyticity of the 
amplitudes in the cut s-plane then leads to |15j 



Nafiis) 



1 



ds' 



Imtaf}[s')Da0{s') 

^{s' -sa){s'-sp){s' - s) 



(19) 



and 



Dap{s) = 1 



[s - So) 



TT 

s - So 



' , Nap{s')Rap(s') 
' (s'-s)(s'-so) 



(20) 



where 



Ka/3(s)r 



Vs - Sg^/s - sp 



•y—TTjK 

(21) 



We have chosen to normalize Nap and Dap such that 
Dap{so) = 1 (a convenient choice that will be employed 
later is sq = 0). The last term in the dispersion relation 
for Dap reflects the so called CDD ambiguity TO]; the 
unitarity relation in Eq.( 14 ) does not uniquely determine 
Dap if iaa vauishcs at some s = Sp^ap, P ~ 1, ■ ■ • , ^p- 
These zeros are then incorporated as poles in D^p with 
7p.a/3 being their residues. It is clear from Eq.(ll) 



that these poles can exist only in the elastic region of 
sk > s > Stt or in the inelastic region s > sk if in- 
elasticity happens to vanish, rj — \ (including the point 
at infinity). At every CDD pole the phase of the elas- 
tic amplitude passes through 180° or the inelastic am- 
plitude vanishes. If the residue of a CDD pole is small 
then Daa{s) will develop a zero on the unphysical sheet 
near the position of the pole, i.e. produce a resonance. 
Thus, in the past it has been proposed to identify CDD 
poles with the elementary quark bound states that turn 
into physical resonances when coupled to the contin- 
uum channels. Indeed it has been shown that in po- 
tential models describing, for example, the scattering of 
a static source with internal structure, the CDD poles 
correspond to excitations of the target [T7]. Asymptot- 
ically, at large s, tap{s oo + ie) < 0{1), and since 



Dap{s ^ oo) = 0(1) it follows from Eqs.([T9| and ([20|, 
that (for P-wave) Napis —?■ oo) = 0(l/sJT~^The set of 



coupled integral equations, Eqs.(|l9|) and (20), gives the 



scattering amplitudes ta,p{s) for all complex s in terms of 
the discontinuity of the scattering amplitudes on the left 
cut and the location of the zeros in the physical region 
(the CDD poles). 

The left hand cut discontinuity plays the role of the 
driving term, which is analogous to the potential in 
nonrelativistic Shrodinger theory and in general it is 
not known. Fortunately, as is clear from Eq.(15), both 



Nap{s) and the production vectors Pq-(s) are real and 
have no singularities in the physical region. Thus it is the 
behavior of the Dap{s) which determine the phase and 
any rapid variation of the isobar amplitudes Fa{s). We 
will use Eqs.( 19 ) and ( 20 1, not as integral equations for N 
and D, but instead we will use what is known about the 
scattering amplitude at the boundary of the right hand 
cut, tapis+ie), with a model for the left hand cut as input 
to determine the denominator functions. Then Eq.(20) 
can be written as an integral equation for D alone 

Dap{s) = 1-n^jp 



(s'-s)(s'-so) 



{s - so) ^^,Dap{s')e-'^^^'^^') sinc^apis') 



(22) 



where (jjap is the phase of tap = \tap{s)\eyi^{i4>ap{s)), 
which has an analytical solution given by 



Dgp{s)^n% 



S Sp^^ct^ 



\Sq - Sg^ap J 

(23) 

The first (second) factor gives the contribution from the 
CDD poles (zeros) and fl is the Omnes-Muskhelishvili 
function, 



^ap{s) = exp - 



s - So 



ds' 



(t^aJ3{s') 



(s'-s)(s'-so) 



(24) 

Phase shifts Sa are determined up to an integer multiple 
of TT and the phase of the amplitude (j)ap is determined 
modulo 2TT. It is customary to remove this ambiguity by 
setting all phase shifts to zero at elastic thresholds, i.e. 
(5Q,(4m^) = 0. This condition is at the origin of zeros 
of Dap being explicit in Eq.(23). With (fyapi^m^) = 
and the asymptotic behavior, Dap{s oo) = 0(1), the 
number of zeros, Nq, and CDD poles, Np, are related by 



cj)ap{00)^TT{Np- Nq). 



(25) 



C. Analytical model for the P-wave amplitude 

If the left hand cut discontinuity of tap{s) were known, 
then the whole amplitude could be reconstructed using 
the N/D method discussed above and the production 
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FIG. 2: Phase shift (top) and inelasticity (bottom) of the 
P-wave nn ampUtude. Data is taken from 19_ (circles) . [20) 
(triangles) , and [21] (squares). The solid line is the result of 
the fit to Stt and r] with the analytical A'-matrix representation 
described in the text. The dashed line is the result of the 
extended parametrization described in Section [II D[ 



vectors Pa{s) could be computed from Eq.(16). Unfor- 



tunately, to the best of our knowledge, only in the case of 
i-^TT is the left hand cut fairly well known [T5]. Thus, one 
needs a model to incorporate the contribution from the 
KK channel. One might as well then construct a model 
that leads to a simple solution of the integral equation 
in Eq.(16|. This is indeed the case if one uses the an- 
alytical if-matrix representation with the typical choice 
of the if-matrix parametrized in terms of simple poles. 
Then the singularity of the scattering amplitude for s < 
is also given by poles and this in turn allows one to solve 
Eq.(16l by algebraic methods. We fix the parameters of 



the 2x2 iiT-matrix so as to reproduce the P-wave tttt data 
from [I5H^ (Fig.[2|; S^^ and r/ are input parameters, and 
the model will give a prediction for 6k- The if -matrix 



parametrization was already used by Haymes at el. to 
interpret their data from [19j . Unfortunately, instead of 



using Eq.(14), the unitarity condition employed in p 
was 

Imiafi{s + ie) = y/s ^ i*^^{s)py{s)0{s - s^)t^/3(s). 

7— 7r,/C 



(26) 
(27) 

(28) 

In contrast, the correct unitarity relation in Eq.([14j) gives 
Im[i-\s)U{s) = -{s - ■s,,)^{l-^-^y,,fj, (29) 



This implies 

Im[i^'^{s)]ap{s) = -(s - 3a)Vs - S^Sap, 

and the if-matrix representation becomes 



which leads to 
where 

■ ^ 1 



/a(s)=/c(0)-- 



s' (s' — s)s' 



(30) 



(31) 



A convenient choice for the subtraction constant, /q(0), 
is to take ReIa{Mp) = 0. Then one of the poles of ifjrTr 
corresponds to the Breit-Wigner mass squared, = 
(0.77 GeV)^, of the p meson. Using the general two-pole 
parametrization of the K matrix. 



4- K 

TTTT — T79 1 T TtTTT, K K 

Alp - S S2- s 



S2 



S2 - S 



IKK 



(32) 



where al = TpMl/{Ml - s^f/'^. By fitting the P-wave 
TTTT phase shift, 5^^^ and the inelasticity, rj, we find T p = 
0.140 GeV, and 

^ = 1.4708 GeV, = 0.199, I3k = 0.899, 

j^K = 0.104, jKK = 1.525, 

(33) 



77r7r = 5.62 X 10 



with the 7's in units of GeV~^. The comparison of 
the phase shift and the inelasticity obtained with this 
parametrization with the data is shown in Figj2] 

Since the K matrix representation of Eq. (|30[) satisfies 
all of the properties of the scattering amplitude discussed 
in Sec II B it is possible to write tajs in the "N/D" rep- 



resentation. We find, choosing to normalize Dafj[s) at 
So = 0, 
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s 



N^^{s) = Kt,-, TT T, D^j,{s) = exp ds' ,/ , ^ n 

[S - SL,l)[S - SLa) \ Js^ S'{s' - S) J 

N^k{s) D^k{s) = 7 r exp / ds r , , 

Nkk{s) = \kk-, r, Dkk{s) =eyi\)[-- I ds ^ 



(34) 



with ^ 5.649, Xkk = 2.271 and X^k = 3.048 GeV^ 
Indeed, as discussed above, the left hand cut is re- 
duced to two poles at sl,i = — 13.87GeV^ and sl,2 — 
—0.787 GeV^, respectively. There are also first or- 
der zeros in N^/^ at Zt^tt = —0.867 GeV^ and zkk — 

— 13.78GeV^. The numerator functions for the elastic 
amplitudes tttt and KK are 0(l/s), and for the inelas- 
tic amplitudes they are super-convergent, i.e. 0(l/s^). 
Asymptotically, as shown in Fig. [3] 07r7r(s — > oo) = 
0(l/log(s)) and (5^ stays below 180°, so there is no 
CDD pole in the tttt channel, which is consistent with 
the Levinson theorem (c/. Eq.(25)). The same is true 
for the KK channel. Above the KK threshold the 
phase of the inelastic amplitude (j^T^K is given by (/Jt^k — 

+ 6k and from the K matrix we find that asymptot- 
ically 0^if(oo) = 27r, which results in two CDD poles 

- one at the p mass, si^ttK = M"^, and the other at 
S2,.K = S2 + P.PkH^k = 3.884 GeV^ 

Having an analytical representation for the scattering 
amplitude enables one to identify the resonance content 
by studying the singularities of t,r7r(s) for s continued 
through the unitarity cuts away from the physical sheet. 
If we define the unphysical sheet II as the one obtained by 
continuing s from above (crossing) the cut < s < sk, 
and sheet III for s continued through the s > sk cut, then 
we find four poles whose location is given in Table [ij The 
p pole is clearly seen as well as the excited p' resonance 
at 1600 MeV that couples primarily to the KK channel. 
The pole on sheet III at 1.1409 - zl).1675 GeV is most 
sensitive to the inelasticity of the KK channel. If we turn 
off the KK channel this pole goes to infinity while the 
positions of the other two remain relatively unchanged. 



TABLE I: Physical poles {^/s in GeV) on sheets II and III. 



II 


III 


0.7638 - i0.0747 


0.7632 - i0.0745 




1.1409-10.1675 




1.6306 - i0.0844 



D. Problems with the /iT-matrix parametrization 

While the K matrix parametrization faithfully re- 
produces the TTTT phase shift and inelasticity data from 
TTTT threshold up to 1.9 GeV, extrapolation beyond this 
range is problematic. The rapid decrease of (pT^T^ around 
s ^ 6 GeV^ seems unphysical and results in an absence 
of the CDD pole at infinity, i.e. ^7r7r(oo) — >■ instead 
of 0^t(oo) — )■ TT [55]. The CDD pole at infinity in the 
elastic tttt amplitude is expected based on the asymp- 
totic pQCD prediction for the pion electromagnetic form 
factor 23J. In the tttt — >■ KK channel, the two CDD 
poles at rrPp and S2 P-kPk Ii-kK a-re clearly an artifact of 
the pole parametrization of the i^T-matrix. A CDD pole 
in the inelastic channel above threshold {e.g the pole at 
S2,-kK — 3.884 GeV^) leads to a discontinuity in a phase 
shift and is unphysical. A pole between tttt and KK 
thresholds is admissible, e.g. the pole at si.ttK = rn?p, 
but its strict overlap with the p mass is an artifact of 
the parametrization. Since the phase space available in 
J/ Tp decay extends up to s^^ ~ 9GeV^ we need to remove 
these unphysical features of the if-matrix amplitude. We 
proceed as follows. The new tttt -> tttt and KK — >■ tttt 
amplitudes will be denoted by i^^™(s) and t^'^'is), re- 
spectively. In the case of the tttt — >■ tttt elastic amplitude, 
we assume that it has a single CDD pole at infinity. We 
thus introduce an effective phase shift and inelasticity 
that asymptotically approach tt and 1, respectively: 



TT + (5„(sk) - 7r)^,s > SK 



(35) 



''^^''^'^ = {l + ivA^K)-in,S>SK 

with ^/sK = 1.9 GeV and 6,^ and 77^ obtained from the 
K-matrix fit below 1.9 GeV (c/. Fig. [5]). The denomina- 
of the effective amplitude 



tor d;'™ 



is then obtained from Eq.(24) with Nq — 
phase, (f>%-fj , given by (see Fig. [3]) 

eff^2iS^ff _ y 

~2i^ 



Nr, 



inew 



= Im In 



n 



(37) 



and 



(38) 
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FIG. 4: Real (upper) and imaginary (lower) part of t^^. The 
dashed line corresponds to the Ji'-niatrix solution of E.q |34| 
and the solid line is the modified _K'-matrix solution, t^™ , 
discussed in Sec. Ill Dl Eq.(|37l). 



FIG. 3: Phase of the tttv (upper) and ttK (lower) amplitude. 
The dashed line is the result of the K matrix parametrization 
from Eq.(34l. The solid line is from the modified K matrix 
parametrization discussed in Sec. |II D| 



For the numerator function A^^^"', we use a simple pole 
approximation to the left hand cut (s^ < 0) 



\new 
s- Sl' 



(39) 



In order to remove the unphysical CDD pole from the 
KK TTTT amplitude for 0^1^(8) in 



for in Eq.([24]), we use (see Fig.[3| 



inew 



2n- 



(40) 



(41) 



with -yi/f = 1.65 GeV. In this case we use the K- 
matrix fit up to a lower energy of 1.65 GeV to be less 



sensitive to the unwanted CDD pole in the K matrix at 
y/S2^TTK — 1.97. There is no effect of this pole in the 
elastic amplitude, and thus for that case we could use 
the ii'-matrix parametrization all the way up to 1.9 GeV 
where data exists. Assuming further that D!^™{s) has 
the same asymptotic behavior as I?^'^™ we add a single 
CDD pole at s""^ m place of the pole at rnj, between 



the TTTT and KK thresholds, 
function we use 



Finally, for the numerator 



Aj-new 



\ new 
S - SL ' 



(42) 



i. e we use the same pole to represent the left hand cut as 
in N^^"". The four parameters A^J^'", A^^J^, sl, and s^'^^j^ 

are determined by simultaneously fitting f"^*" and t"^™ to 
TT and K phase shifts and inelasticity in the range 2mj^ < 
y/s < 1.9 GeV and 2mK < \/s < 1.65 GeV, respectively. 
The comparison with the /f-matrix solution is shown in 
Figs. |4] and [5] and the fit yields A^J^"" = 0.750, A^fff = 



0.0477, Sl 



-1.328 GeV^ and s 



1,7tK 



0.220 GeV^ 



As expected, the location of the left hand side pole falls 



4 6 
s [GeVI 
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FIG. 5: Same as Fig. ^ for Ui 
iZ'i^ from Eq.ra. 



from Eq. ( 34 1 (dashed) and 



FIG. 6: Real (dashed) and imaginary (solid) part of the 
inverse of D"^™(s) (upper) and D"'j^{s) (lower) used in the 
computation of the isobar form factor, cf. Eq.(44|. 



between the two left hand side poles of the iC-matrix 
parametrization. In Fig. [6] we show the inverse of the 
denominator functions Z?^'^*^' and D^'j^ . 



E. Interpretation of the J/i/i — > Stt data 

With the left hand cut singularities of the scattering 
amplitudes given by a simple pole, {cf. Eqs.(39),(42)) 
from Eq.([l7| it follows that ImPa{s) = 0. Thus Pa(s) is 
analytica in the entire s-plane and therefore given by a 
polynomial, 



P„(s) = (s-Si)Ca(s). 



(43) 



The first term is responsible for removing the left hand 
cut singularities from N^™{s) and making Fa{s) in 
Eq.(IT5| analytical for s < 0. The bound \Pa{oo)\ < 1 



restricts Cq(s) to be at most a first order polynomial in 
s. Thus the final solution to Eq.(15l has the form 



1 + 



1 + gks 



(44) 



The first term corresponds to J/ip — > (7r7r)p7r and the 
second to the re-scattering contribution from J/i/j — 
{KK)pTi — > (7r7r)p7r. 

The Dalitz distribution of Btt events from J decays is 
shown in Fig. [8] and the striking feature is the depletion of 
events in the center of the plot. This is to be compared 
with the distribution shown in Fig. |9j which has been 
generated with r.^K = 0. The three bands originate from 
the p meson contribution to and the large contri- 

bution from the p'(1600) resonance leads to a significant 
population of events in the middle of the Dalitz plot that 
is not seen in the data in Fig. [8) Furthermore in the 
data there is a large contribution near the tails of the p 
bands, which are absent if only the direct Btt production 
is considered. We thus consider the full amplitude from 



Eq.(44) and float the three parameters a.,^,aK and Tt^k 
to obtain a distribution that best resembles the data. We 
find little sensitivity to the term proportional to aK and 
thus set aK = 0. The parameter is relevant since it 
controls the tail of the p resonance and so is r^^K which 
determines the relative strength of the KK contribution 
which interferes with the tttt amplitude in the p'(1600) 
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using the same parameters as in Fig. |10| S GeV^ 




FIG. 9: Dalitz plot distribution with the single nn channel 
only i.e. -Fi(s) = qTrPTr/Dn instead of Eq.(44|. 



FIG. 8: The J/i/) -5> 7r+7r"7r° Dalitz plot distribution from 
the BBS Collaboration [2]. 




3 4 5 
S GeV^ 



region and reduces the contribution at the center of the 
Dalitz plot. In Fig. 10 we show the event distribution 
-1.5 X lO^GeV"^ and r„/K = -1.3 x 10"^ 



using 



The normalization constant N is at this stage arbitrary 
since we are not determining the absolute value of the 
branching ratio. 

Now, inspecting the Dalitz plot in FigJTO] and the 
plot of the function in Figj?} it is seen that the 

KK channel can indeed bring theory closer to the data 
by enhancing the tttt contribution in the energy range 
1 GeV < ^/s < 1.5 GeV and reducing the strength of the 
p'(1600) peak. 



FIG. 10: Dalitz plot distribution from Fi in Eq.(44l with 
both the TTTT and KK channels with — —1.5 x 10~ GeV~^ 



and r. 



tt/K 



-1.3 X 10" 



III. SUMMARY 

We have studied the effects of inelastic tttt scattering 
on the J/ip — T' Stt Dalitz plot. We have seen that the 
KK — > TTTT channel can significantly alter the shape of the 
Dalitz plot, especially at higher tttt masses. This brings 
the observed data closer to the phenomenological expec- 
tations based on tttt P-wave scattering. These coupled 
channel effects will become even more important as ex- 
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perimental data sets grow larger, for example at BES III, 
where 1 billion J/?/' decays are expected. 
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